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The light-cone distribution amphtudes for the axial current are derived within the instanton vac- 
uum model (IVM), which incorporates nonperturbative effects including spontaneous chiral sym- 
metry breaking. This allows to extend applicability of the dipole approach, usually used in the 
perturbative domain, down to Q^ —>■ 0, where partially conserved axial current (PCAC) imposes a 
relation between the neutrino-production cross section and the one induced by pions. A dramatic 
breakdown of the Adler relation ( AR) for diffractive neutrino-production of pions, caused by absorp- 
tive corrections, was revealed recently in pj. Indeed, comparing with the cross section predicted by 
the dipole phenomenology at Q^ — >■ on a proton target we confirmed the sizable deviation from the 
value given by the AR, as was estimated in [l[ within a simplified two-channel model. The dipole 
approach also confirms that in the black-disc limit, where the absorptive corrections maximize, the 
diffractive cross section ceases, on the contrary to the expectation based on PCAC. 
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I. INTRODUCTION 



Due to its V-A shape the neutrino-hadron interactions 
possess a very rich structure. However, because of the 
smallness of the cross-sections until recently the experi- 
mental data have been scarce, mostly being limited to the 
total cross-sections. With the launch of the new high- 
statistics experiments like MINERi^A at Fermilab [2], now 
the neutrino-hadron interactions may be studied with a 
better precision. The V-A structure of the neutrino-quark 
vertices enables us to study simultaneously (VV), (AA) 
and (VA) correlators in the same process. 

The properties of the vector current have been well stud- 
ied in the processes of deep inelastic scattering (DIS) of 
leptons on protons and nuclei, deeply virtual Compton scat- 
tering (DVCS), real Compton scattering (RCS), and vec- 
tor meson production. The standard approach in the de- 
scription of such processes is based on the large-Q^ factor- 
ization of the cross section into a process-dependent hard 
part, which is evaluated in pQCD, and a universal target- 
dependent soft part. The latter is extracted from fits to 
experimental data. Factorization, however, is not valid at 
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small photon virtualities, where one can rely on the dis- 
persion relation, or on the assumption of generalized vec- 
tor meson dominance (GVMD) 3 6]. Such a description, 
however, involves a lot of ad-hoc modeling. 

An alternative phenomenology for high-energy QCD pro- 
cesses is based on the color dipole approach [3]. One as- 
sumes that before interaction, the projectile (virtual W- or 
Z-boson in case of the neutrino scattering) fluctuates into 
a quark-antiquark dipole. After the dipole is formed, it 
scatters in the field of the target and then fluctuates back 
to the final hadron Q. Recently the color dipole approach 
has been successfully applied to the description of different 
reactions with vector currents (see |8l-l27| and references 
therein) . 

For the axial current the situation is more complicated, 
especially at small Q^, because the chiral symmetry break- 
ing generates the near-massless pseudo-goldstone mesons 
(pions) . Straightforward extension of the vector dominance 
model to the axial current leads to the so called Piketty- 
Stodolsky paradox [28, ^], which appears because axial 
meson dominance is broken by a large contribution of multi- 
pion singularities in the dispersion relation [29l . |30| • The 
dipole description is free of this problem, because in this 
model there is no explicit hadronic degrees of freedom and 
the interaction occurs via dipole scattering. 

According to the Adler theorem [3l|, JS^], based on the 
hypothesis of partial conservation of axial current (PCAC), 
the neutrino-proton interactions cross-section at zero Q^ 



is proportional to the cross-section of the pure hadronic 
process, where the heavy intermediate boson is replaced by 
a pion, 
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where F denotes the final hadronic state; y — v/E^] Ey is 
the energy of the neutrino; and v is the energy of the heavy 
boson W , or Z , in the target rest frame. The chiral sym- 
metry is vital and should be embedded into any dynamical 
model which is used for calculation of the cross section at 
small Q"^. 

In what follows we entirely neglect the lepton mass (accu- 
rate for neutral currents or electrons), which can be easily 
incorporated Q , but is dropped for the sake of simplicity. 
If one interpreted the Adler relation (AR) (Eq. [Ij in terms 
of pion pole dominance, one would arrive at a vanishing 
cross section. Indeed, the pion pole term in the amplitude 
contains a factor g^ , which multiplied by the conserved lep- 
ton current /^ terminates this contribution [l|, |28| . 123 . |33| . 
Other, heavier meson states provide a final contribution, 
but have to conspire to mimic the pion pole term, as is 
dictated by the PC AC relation Eq. ([T]). 

Such a fine tuning looks miraculous if one has no clue of 
the underlying dynamics. Similar paradox is known for the 
1/Q^ behavior of DIS cross section. In QCD this is known 
as a result of color screening [3] , leading to the effect called 
color transparency. While this is rather obvious within the 
color dipole description, it becomes extremely sophisticated 
in hadronic representation. Indeed, expanding the current 
over hadronic states one may arrive at a problem called 
Bjorken paradox 34]. Namely, how it happens that many 
hadronic states, having large sizes and large cross sections 
conspire in a way that all together they act like a small 
hadron with a tiny, ^ 1/Q^ cross section. The solution is 
known, the off diagonal diffractive amplitudes are negative 
and cancel with diagonal ones [3a|- This however cannot 
be proven, unless one employs an explicit model describ- 
ing the features of the hadronic states and the diffractive 
amplitudes. This is why the color transparency effect has 
not been understood within the GVMD, but was revealed 
in the color dipole representation [7] . 

Similarly, in order to test the mysterious relation between 
the contribution of heavy hadronic fluctuations and pion, 
one should switch to the dipole representation and employ 
models for the distribution amplitudes (DA) of the axial 
current which have built-in chiral symmetry. Recently, we 
used the DA of the vector current calculated in the IVM for 



the evaluation of several processes |ll| - KL4| . In this paper we 
apply the IVM to construct the DAs for the axial current 
and pion and use them to calculate the neutrino-production 
cross sections. Since the IVM includes spontaneous chiral 
symmetry breaking, the qq DAs of axial current and pion 
should automatically satisfy PCAC, and in the small-Q^ 
limit reproduce the Adler relation ([IJ. 

Notice that the color dipole description is valid only at 
high-energies or at small xb ^ 1, where the contribution 
of quark exchange (reggeons) is suppressed as l/y^- For 
moderate energies another mechanisms, such as e.g. for- 
mation of resonances in the direct channel [36l |37[ , and/or 
reggeon exchange in the crossed channel may be impor- 
tant [23, [33 • In this paper we do not consider those correc- 
tions, but concentrate on the well developed small-x dipole 
phenomenology. 

Experimentally, the neutrino-production of hadrons on 
protons and nuclei has been studied in the recent exper- 
iments K2K '4^-42], MiniBoonepl SI, NuTeV [H^i^ 
(see also review 29] and [30, l48l - [5l| for references to earlier 
neutrino experiments). For high energy neutrino scatter- 
ing, there are data from the early bubble chamber experi- 
ments [33.1^] with energies up to 100 GeV, though with low 
statistics and only for the total (integrated) cross-sections. 
Currently, with the launch of the high statistics experiment 
Minerva at Fermilab [3, [23| , the precision of measurements 
should be considerably improved, and data for the differ- 
ential cross-sections at high energies will become available. 

In this paper we consider a particular process - diffrac- 
tive single pion production on a proton target. As was 
demonstrated in [l| and below, this process provides a most 
sensitive way to test PCAC in high energy neutrino inter- 
actions. Besides, it generates an important background to 
the measurements of neutrino oscillations [53|-|56j, and is 
also important for the neutrino astronomy of astrophysical 
and cosmological sources. 

The paper is organized as follows. In Section |TT] we 
present the color dipole formalism. In Section IIIII we per- 
form calculations of the DAs of the axial current and pion. 
In Section IIVI wc calculate the overlap of the DAs for the 
axial current and of the pion and found it to be proportional 
to 9^1, what terminates this contribution to the neutrino- 
production of pions due to conservation of lepton current. 
In Section |V] we present the numerical results and summa- 
rize the observations in Section IVll 



II. DIFFRACTIVE PRODUCTION OF PIONS 

The cross section of diffractive neutrino-production of a 
pion on a proton, vp — >■ lirp^ has the form, 
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where win is the nucleon mass; L^^ is the lepton tensor; 
and W^ is the amphtude of pion production by the axial 
current on the proton target. In the color dipole model this 
amplitude has the form 






dMMrid'r2^"{P2,r2) (3) 



X ^'^(/3i,ri;/32,r2;A)*;^(/3i,ri) 

where ^'^ and, ^1^ are the distribution amplitudes (DAs) 
of the pion and axial current respectively, and A'^{...) is 
the dipole scattering amplitude. The axial current DA '^^ 
contains a pion pole, whose contribution to the amplitude 
is proportional to g^, because the pion is spinless. This 
factor terminates the pion pole because of conservation of 



the lepton current. As we assumed, the lepton is massless, 
otherwise the pion pole contribution is not zero and leads 

to corrections of the order of O 






The amplitude A'^ {(3i,ri; I32,r2', A) in ^ depends on the 
initial and final quark transverse separations ri.2, fractional 
light-cone momenta /3i^2, and transverse momentum trans- 
fer A. This is a universal function dependent only on the 
target but not on the initial and final states. In addition 
to the axial current contribution, in ([3]) there should be 
the contribution of the vector current. This contribution 
involves a poorly known helicity flip dipole amplitude Ad, 
which vanishes at high energies as l/v. Besides, at small 
Q^, the vector current contribution is suppressed by a fac- 
tor Q^. At high energies, in the small angle approximation, 
A/y^ <C 1, and the quark separation and fractional mo- 
menta /3 are preserved, so 

A''{l3i,ri]l32,f2;Q^x,A) « 5(/3i - /Sa) <5(ri - r2) (4) 

X (e + *)3m/g^(f,A,/3,x) 

where e is the ratio of the real to imaginary parts, and 
for the imaginary part of the elastic digole amplitude we 
employ the model developed in [111 . |57h59| , 



J 



'^m f^q [r, 



N 1^ ^ 



p,x] = — - — exp 



B{x) Rljx) 
2 16 



A 



r 



,-zfir-A _^_ g»(l-,9)r.A __ 2e<5 -/3)*='-Aj 



(5) 



The phenomenological functions cro{x), Rq{x) and B{x) are 
fitted to DIS and p electroproduction data. We rely here on 
the Bjorken variable x — Q^ /2(pq), which has the meaning 
of fractional light-cone momentum of the parton only at 
large Q^. At low Q^ important for the axial current, one 
should switch to an energy dependent parametrization, as 
is explained in Section FVl 



For the forward scattering, A — ^ 0, the imaginary part of 
the amplitude (|4]) reduces to the saturated parameteriza- 
tion of the dipole cross-section proposed by Golec-Biernat 



and Wiisthoff (GBW) 8], 



A = 0, /3, a; 
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1 — exp 
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(6) 



Generally speaking, the amplitude /i^(...) involves nonper- 
turbative physics, but its asymptotic behavior at small r is 
controlled by pQCD [7|: 

up to slowly varying factors ^ ln(r) [7|. 

Calculation of the differential cross section also involves 
the real part of scattering amplitude, whose relation to the 



imaginary part is quite straightforward. According to [60[, 
if lim ( ™-' ) is finite, then tlie real and imaginary parts 

s— >-oo V * / 

of the forward ampHtude are related as 



7^e/(A==0) = s"tan 



91ns 



3m /(A = 0) 

s^ ' 

(7) 

In the model under consideration the imaginary part of the 
forward dipole amplitude indeed has a power dependence 
on energy, Xto/(A — 0; s) ^ s"^~^, where ajp is the in- 
tercept of the effective Pomeron trajectory. Then Eq. ([7]) 
simplifies to 

1^ =tan(f(a^-l))^e. (8) 

This fixes the phase of the forward scattering amplitude, 
which we retain for nonzero momentum transfer, assuming 
similar A-dependences for the real and imaginary parts. 



III. DISTRIBUTION AMPLITUDES AND THE 
INSTANTON VACUUM MODEL 

In this section we define the DAs and give a brief de- 
scription of the instanton model used for their evalua- 
tion (see |6ll - [63j and references therein). 



A. Instanton vacuum model 

The central object of the model is the effective action for 
light quarks in the instanton vacuum, which in the leading 
order in Nc has the form 162, e3I 



S = f d'^x(ylnX + 2^^{x) (9) 

- 'ip{p + v + aj5~m- cLf » $ • T™ » fL) ip) 

where Tm is one of the matrices, Tm — 1, ir, 75, JT75; ip and 
$ are the fields of constituent quarks and mesons respec- 
tively; N/V is the density of the instanton gas; to w 5 MeV 
is the current quark mass; v = v^,^^ is the external vector 
current corresponding to the photon. L is the gauge factor 
defined as, 



L{x,z) = Pexp ( z I dC^K(C)+ap(C)75) ] • (10) 
L{x,z) = 7oL^(x,z)7o (11) 



It provides the gauge invariance of the action, and f{p) is 
the Fourier transform of the zero-mode profile in the single- 
instanton background. In this paper we used for evalua- 
tions the dipole-type parameterization [621 



f{p) 



L' 



L^ 



pz 



with L- 850 MeV. 

In the leading order in N^^ we have the same Feynman 
rules as in the perturbative theory, but with a momentum- 
dependent quark mass /i(p) in the quark propagator 



S{P) = 



1 



(12) 



p — j-i{p) + iO 

The mass of the constituent quark has the form 

H{p)^m + Mf{p), 

where to, w 5 MeV is the current quark mass, M w 350 MeV 
is the dynamical mass generated by the interaction with 
the instanton vacuum background. Due to the presence of 
the instantons the vector current - quark coupling is also 
modified. 



vt,r 



-^V = v- 
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In addition to the vertices of the perturbative QCD, the 
model contains the nonlocal terms with higher-order cou- 
plings of currents to mesons. The exact expressions for the 
nonlocal terms \/"°"' ^ Anoni (jgpejK;^ on the choice of the 
path i n (1101). so one can find different results in the liter- 
ature [6J-l67|. However, for the longitudinal parts of the 
axial and vector currents important here, this ambiguity 
cancels out and couplings have the form 
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where pi, p2 are the momenta of the initial and final quarks. 



B. Axial current distribution amplitudes 

The distribution amplitudes of the axial current are de- 
fined via 3-point correlators 

^/3 ^ y d^i €-'"< (0 1^ (y) rV' {x) JliO I 0) , (13) 



where Ja{^) is the axial isovector current and F is one of 
the Dirac matrices. Due to the spontaneous chiral symme- 
try breaking and existence of the near-massless pions the 
hadronic structure of the axial current differs from that of 
the vector current. In particular, the axial current may fluc- 
tuate into the pion state before production of the qq pair. 
Thus the correlator (|13p has two contributions, schemati- 
cally shown in the Figure [TJ 







FIG. 1: The distribution amplitude has two contributions, with 
intermediate heavy axial states and with a pion, labeled with 
(bulk) and (pion) respectively. 

One term corresponds to the combined contribution of 



the intermediate heavy states (oi meson, Stt, etc.), another 
one corresponds to the axial current fluctuating into a pion. 
Due to the built-in chiral symmetry, the two contributions 
are connected by PCAC, so for the full DA we have 



This form of the DA reflects the relation between the pion 
pole and the bulk of heavy states contribution imposed by 
PCAC, which has been discussed above. In what follows we 
concentrate on the part of the amplitude presented in the 
dispersion relation by the bulk of heavy states excluding 
the pion pole (left pane in the Figure [1]), tacitly assuming 
that the full distribution amplitudes are determined using 

m- 

For the part of the axial current presented by the bulk 
of heavy states, the DAs may be defined similar to the 
distribution amplitudes of the axial meson [68[ : 
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(17) 



(0|Vi(y)75^(x)|A(q)) = /le(^'-W d/Je't'^P-^+^P-)^^, 



r 



(18) 



where q is the momentum carried by the axial current; /3 is the fractional light-cone momentum; /3 — 1 — /3, 



e(-'') ^ e^'^''(g) is the polarization vector of the axial me- 
son with polarization state A; z — x — y; p^ is the "positive 
direction" vector on the light-cone; n^ is the "negative di- 
rection" vector on the light-cone. Light cone vectors p, n 
are chosen in such a way that the vector q does not have 
transverse components. The normalization constant /^ is a 
dimensional parameter introduced in order to make the dis- 
tribution amplitudes dimensionless. Its value is fixed from 
the condition 



dp^jp) = 1. 



We defined an "effective" axial state [^(^^(q)) as 



(19) 



A^^\q)) = I d^xe-^^--ef\q)Jl{x) |0) . (20) 



The DAs have the following twists: $|| (^), $^(/3) are twist- 

2; g±^ ) ffj^ ' '''I I J ^M ^^^ °^ twist-3; 53, /13 are of twist-4. All 
the wave functions in (|15ll6p are chiral even; all the wave 
functions in (|17ll8p are chiral odd. 



C. Pion distribution amplitudes 

A sp inless pion has only four independent DAs defined 
as [si]: 

(0|V5(2/)7m75^(x)|^(9)) = lU t d/3e^(^^-^+^P-^) (21) 
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Twist counting is the following: 4'2;tt is a single twist-2 
function (it was evaluated earlier in |70l - [7a |). (pf.l and (fy^.l 
are twist-3, ■04;7r is the twist-4 DA. 



The fuU expressions for the DAs ([I5l)-(IIH1) and (^Tll- ip^ 
are given in Appendixes lA II and IA2I respectively. 

The DAs for the vector current are presented in Ap- 
pendix IA3I However, the vector current does not con- 
tribute to the color dipole amplitudes of pion production. 
Although it contains nonzero components, their overlap 
with the pion DAs is zero. The vector part vanishes because 
the color dipole amplitude does not fiip helicity. Within the 
VDM approximation such components may be expressed 
via the pN — > iiN scattering amplitudes, which exist only 
due to quark-antiquark (Reggeon) exchange in the cross 
channel. This is beyond the employed dipole phenomenol- 
ogy corresponding to gluonic (Pomeron) exchanges. 



IV. DISAPPEARANCE OF THE PION POLE IN 
THE DIPOLE REPRESENTATION 

As was emphasized above, the pion pole contribution to 
the pion production amplitude vanishes because of lepton 
current conservation (up to the lepton mass). This nontriv- 
ial observation of [l|, [2^ [2^, jSj] is in variance with the naive 
interpretation of the AR Eq. ((T|), which relates diffractive 
neutrino-production of pions, v+p — > ? + 7r+p, with elastic 
pion-proton scattering, it -\- p ^ tt -\- p. It is tempting to 
interpret this relation as pion pole dominance, i.e. neutrino 
fluctuates to a pion, which then interacts elastically with 
the proton target. If this were true, the amplitude should 
be maximized in the so called black disk limit, which cor- 
respond to unitarity saturation when the imaginary part 
of the partial elastic amplitude reaches the maximal value 
allowed by the unitarity relation. 

On the other hand, if the pion pole does not contribute 
[H, [23, [2E IsSl a^s is stressed above, all hadronic fluctua- 
tions of the neutrino contributing to vl/(^""=) are heavier 
than a pion, so all diffractive hadronic amplitudes of pion 
production are off-diagonal. Such amplitudes vanish in the 
black-disc limit, so the pion cannot be produced diffrac- 
tively. The source of such a dramatic breakdown of PCAC 
was identified in Jj] as a result of strong absorptive correc- 
tions. Of course the deviation from the PCAC prediction, 
AR, on a proton or nuclear targets, which may be far from 
the unitarity bound, is not so dramatic, as was calculated 
in ,1]. 

In this section we present an explicit demonstration of 
disappearance of diffractive pion production in the black- 
disc limit relying on the dipole model. Namely in this 
regime all the partial elastic amplitudes ([S]) reach the uni- 



tarity bound, so become independent of the dipole trans- 
verse separation r, and the equation ^ simphfics to just an 
overlap of the initial (axial current) and final (pion) light- 
cone wave function. We intent to demonstrate that this 
overlap vanishes. 

The amplitude of pion production in this regime has the 
form, 
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X ^i^) (/3, r; q - A) vj/^ (^^ ^. g) ^ (34) 
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where the index a numerates all the distribution ampli- 
tudes. This is suppressed, since transition from spin-1 to 
spin-0 requires helicity flip for one of the quarks in the 
quark-antiquark pair. Now we would like to demonstrate 
explicitly that such suppression indeed takes place in case 
of the perturbative QCD model. 

The distribution amplitude of the meson state is defined 
as 



dze'^'e''^''{0 |^(z,r)rV'(0)| M{q)) = $m (/3,r;g) 



where the separation between the quark and antiquark has 
a "minus" and transverse components {z,r), and F is one of 
the Dirac matrices (l, 7^, 7/j, 7,^75, cr^;/) multiplied by the 
proper isospin factor (for isoscalar this is one, for isovector 
mesons it is ir etc.). The exact expression on the right 
hand side depends on the matrix F, spin of the meson and is 
usually given as a twist expansion over all possible Lorentz 
structures which may be constructed from T,q,f and the 
polarization vector e{q). 



imj^ 



FIG. 2: Diagram corresponding to the distribution ampli- 
tude dm in pQCD 



In the leading order of as (which is justified for very large 
q'^) the corresponding DA may be represented as a simple 



(25) diagram shown schematically in Fig. [51 Then we have 
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(26) 



(2/c+fc- - /ci - m2 + iO) (2 (fc+ - q+) (fc- - q-) - {kj_ + q^f - m^ + iO 



where the function /(fc, q) = 

(lk + m]rik + q + Tn] Tm ) depends on the spins 

of mesons and Dirac matrices and is not important for a 
moment. 



Taking the integral over k , we get 

d^kS (/3 - k+/q+) e'^^^'^Tr [S{k)TMS{k - q)T] (27) 
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2{l3{l-P)q^-{kl+m^j) 
Straightforward evaluation of the overlap of the two func- 



tions {^At^tt) is quite tedious, however we may signifi- 
cantly simplify the evaluations using completeness of the 
Dirac matrices, viz. 



E-p(n)-p(n) 
^ ap^ a' 13' — "aa 



Saa'Sa 



(28) 



so the product of numerators of the two DAs is now con- 
verted to the effective diagram shown in the Fig. ([3]). 



^ * — > ^ 



FIG. 3: Diagram corresponding to overlap of the distribution 
amplitudes 



Straightforward evaluation gives 
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Qf, (4to2 + A^) - A^ (4m2 + 4k ■ q + q"^) + 2fc^ (-2fc ■ q + A"^ - q'^ + q ■ A) 



4 (/3(1 - /3)g2 _ {kl + m2)) (/3(1 - /3) {q - A)^ 
I 



(kj 



As we can see, the result is proportional to 0{m) ~ O (m^) 
and thus is suppressed in the chiral limit. We expect that 
the same result is valid for the nonperturbative DAs eval- 
uated in the instanton vacuum model. 



The function ao{s) is fixed by the condition 



dVa,,(r,s) / d/3|vl/- (/3,r)|2 = aZ{s). (32) 



V. NUMERICAL RESULTS 

The Bjorken variable x = Q^/ (2p • q), used at high Q^, 
is not appropriate at small Q^, where it does not have the 
meaning of a fractional quark momentum any more, and 
may be very small even at low energies. For the case Q^ = 
0, where the AR holds, x defined in this way would be 
zero. Therefore, one should rely on the phenomenological 
dipole cross section which depends on energy, rather than 
X. At small Q^ we employ the s-dependent parametrization 
of the dipole cross section [TJ], which is similar to the x- 
dependcnt GBW parametrization \8\, but is more suitable 
for soft processes 



Ro{s) 



1 



Ms) ( 

0.88 /mx (—] 



rVKis) 



0.14 



(30) 
(31) 



These parameters and the scale sq = 1000 GeV^ are fitted 
to data on DIS, real photoproduction and irp scattering. 



A. Corrections to the AR 

It was pointed out in \Y\ that the AR ([T|) applied to 
diffractive neutrino-production of pions, should be broken 
by absorptive corrections, which affect the left- and right- 
hand sides of Eq. ^ differently. Therefore, the AR Eq. ^ 
cannot hold universally, since the magnitude of absorptive 
corrections is target dependent. The corrections reach max- 
imum in the black-disc limit (e.g. on heavy nuclei), where 
the AR is severely broken as was demonstrated above in 
Sect. II VI It was revealed in T] that the AR is not accurate 
even on a proton target, a deviation of about 30% was esti- 
mated for diffractive pion production on a proton. Unfor- 
tunately, the baseline for comparison of the left and right- 
hand sides of Eq. ([T|) is ill defined, it is not clear whether 
the AR should hold with or without absorptive corrections. 

The dipole phenomenology, which is adjusted to data, is 
free of this uncertainty, it does not need to be corrected for 
absorption. One can calculate the left and right-hand sides 
of Eq. ([T]) on the same footing. Important is to use the same 



crJ^'J(s) in ([5^ . as in the right-hand side of ([T]). We employ 
only the pomeron part of the cross section parametrized as 
arjJ(s) = 23.6m6x(s/so)°o«. 

Calculation of the left and right-hand side of Eq. ([T]) 
clearly demonstrates that absorptive corrections affect 
them differently. Indeed, the amplitude in the left-hand 
side of (HI) is given by Eqs. © and ([3]). The first term in the 
dipole cross section (15(11) is independent of r, therefore its 
contribution to neutrino-production amplitudes vanishes, 
as was demonstrated in Sect.|IVl So the production ampli- 
tude is suppressed by the second exponential term in (PU)) . 
This suppression becomes stronger with energy, since R{s) 
decreases, and at very high energies in the Froissart regime 
the left-hand side of ([T]) vanishes. At the same time, the 
pion-proton cross section Eq. p2l) is dominated by the first 
term in (|30p and reaches maximum in the Froissart regime. 

Now we are in a position to evaluate the accuracy of the 
AR for diffractive neutrino-production of pious on protons. 
In Figure U] we plotted the ratio of the cross-sections cal- 
culated with the color dipole model (left-hand side of ([1])) 
and using the AR (right-hand side of ([Ij), 



Kar{s) 
as was defined in [l|. 



da dipole /dtdvdQ'^ 



doAR/dt dv dQ^ 



(33) 
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FIG. 4: Ratio of the cross-sections calculated within the color 
dipole model and using the AR ([!]) at Ax = 0. 

As was expected, Kab. < 1 due to different structures 



of the absorptive corrections to diffractive pion production 
and elastic pion scattering cross-sections. The deviation of 
Kar from unity is significant, even somewhat larger than 
was estimated in [l|. The ratio is falling at high energies 
towards the Froissart limit, where it eventually vanishes 
when Ro{s) — >■ 0. 



B. Predicted cross sections 

Most of the data on neutrino-production of pious on pro- 
tons have been available so far only at energies close to the 
resonance region [2^. Data at higher energies are scarce 
and have rather low statistics [3J, [S^l ■ Because the dipole 
formalism should not be trusted at low energies, we provide 
predictions for the ener gy r ange of the ongoing experiment 
Minerva at Fermilab 0, [23| • 

The Q^ dependence of the diffractive cross section de- 
serves special attention. It would be very steep at small 
Q^, if the pion dominance were real. However, since the 
pion pole is terminated due to conservation of the lepton 
current, the Q^ dependence is controlled by heavier sin- 
gularities. In the approximation of an effective singular- 
ity at Q^ — — M^ [l| one should expect the dipole form 
ex {Q^ + NP)~'^ . Within the dispersion approach the effec- 
tive mass scale M is expected to be of the order of 1 GeV 
[H, m, HI. Within the dipole description the Q^ depen- 
dence is controlled by the the IVM mass scale, which is of 
the order of 700 MeV. 

In Fig. [SJ we plot the forward diffractive neutrino cross 
section scaled by the factor (Q^-l-M^)^, where the parame- 
ter M is adjusted in a way to provide a flat Q^ dependence 
at Q2<2GeV2. 

Indeed, we found that at M = 0.91 GeV the scaled cross 
section is constant up to rather large Q^ ~ 3GeV^, but 
substantially deviates from the dipole form at larger Q^ . 

The t-dependence of the cross section is controlled by the 
employed model Eqs. (0 for impact parameter dependence 
of the dipole amplitude. The results for i-dependence of 
the invariant cross-section are shown in Fig. [B] for several 
fixed values of Ei, and Q^ = A GeV^ . For this calculation 
we fixed y — 0.5. 

The forward invariant cross-section Eq. ([2]) of diffractive 
neutrino-production of pions on protons is depicted in the 
Fig. d?]) as function of v at several fixed values of y and Q^ . 

These calculations performed in the dipole approach are 
controlled in Eq. ([3]) by the light-cone DAs of the axial 
current and pion, which we extended to the soft interaction 
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FIG. 5: The Q^-dependence of the cross section of diffractive 
neutrino-production of pions scaled by factor {Q^ + AI^)^ at 
neutrino energy Ei, = 20 GeV and different y. The mass pa- 
rameter M — 0.91 GeV is adjusted to minimize the variations 
of the scaled cross section at small Q^ . 



regime basing on the instanton vacuum model. It worth 
reminding that neither s-channel resonances, nor reggeons 
are included in the parametrization ([30|) of the universal 
dipole cross section, so the results are trustable only at 
sufficiently high energy v. 

Experimental data for neutrino-production cross section 
are usually presented as function of neutrino energy E^, in- 
tegrated over v. Unfortunately, in this form one cannot 
separate physics of low and high energies. Indeed, the in- 
tegration over v results in the finite contribution of small- z^ 
region, which is dominated by s-channel resonances. This 
small-:/ contribution is constant at any high neutrino en- 
ergy El, and its magnitude is comparable with the diffrac- 
tive part. 

Usually in low statistics experiments one integrates 
the multi-dimensional distributions presenting the re- 
sults as function of one variable. As such a variable 
we chose the cm. energy of the diffraction process, 
Q^ + 2mNV. Then we calculate the W- 



W ^ ^r 
distribution as. 



'AT 



^^2W [ dvdt dQH (W^ -{p + qf 
dW J 



da 



diydtdQ^' ^^^^ 
In the left pane of the Figure|8]we plotted M^-dependence 
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FIG. 6: The t-dependence of the cross section of diffractive 
neutrino-production of pions at different Q^ for neutrino energy 
El, = 20 GeV and y = 0.5. 



of the cross-section p4p for several fixed values of E^, . We 
see that the W^-dependence significantly varies with iJ^, 
therefore one should average the cross section Eq. (p4| 
weighted with a realistic neutrino energy distribution. 



da_\ 



dE„ p {E„) 



da 
dW' 



(35) 



where the neutrino spectrum p (E^) is normalized as 

fdE,p{E,) = l. (36) 

As an example, we performed calculations with the neu- 
trino energy spectrum of the MINERvA experiment [2[. 
We considered three different _Ei,-distributions correspond- 
ing to low (LE), medium (ME) and high energy (HE) beam 
configurations. The results are depicted in the right pane 
of the Figure [51 

We also compared our results for the VF-distribution of 
neutrino diffi'active events with data of the WA21 experi- 
ment at CERN [54I . We performed averaging over neutrino 
energy with the spectrum p{Eu) given in [74|- The results 
are depicted by dashed curve in Fig. [HI Since the low en- 
ergy region is affected by reggeons, which we have neglected 
so far, we added their contribution to the dipole cross sec- 
tion (1221) ■ 
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FIG. 7: [Color online] Forward neutrino-production cross- 
section of pions as function of f at several fixed values of y 
and Q^. 



aZis) = [13.6s°o« + 19.2s-°-4V^ 



(37) 



Q^, where the nonperturbative effects are unavoidable. We 
calculated the distribution amplitudes for the axial current 
(Sect. lllTB|) and for the pion (Sect. llirC|) on the same foot- 
ing, within the instanton vacuum model (Sect. IIII Ap . The 
model possesses the chiral symmetry properties, what guar- 
antees a correct, controlled by PCAC behavior at small Q^. 

Although the dipole approach does not involve explic- 
itly the intermediate hadronic states, absence of the pion 
pole can be tested in the "black-disc" regime, where all the 
partial elastic amplitude saturate at the unitarity bound. 
Indeed, the direct calculation performed in Sect. IIVI con- 
firmed that diffractive pion production ceases, what may 
happen only if the pion pole does not contribute. 

The dipole description also offers an unbiased way to test 
the AR on a proton target. This relation is expected to be 
broken by absorptive corrections [l|, which are implicitly 
included in the phenomenological dipole cross section. We 
rely on the dipole cross section parametrized in the satu- 
ration form Eq. ((30|. well confirmed by data for electro- 
magnetic processes |8|. We found a significant, about 40% 
deviation from the AR on a proton target (see Fig. S]). 

A much stronger breakdown of the AR is expected for 
nuclei [l|, and we plan to evaluate those effects employing 
the techniques developed here. 



The result shown by solid curve describes the data much 
better. For comparison wc plotted also the prediction 
based on the AR with the realistic pion-proton cross section 
Eq. 1321). 



VI. SUMMARY 

We developed the dipole description for high-energy neu- 
trino interaction, in particular at low Q^, whether the 
PCAC hypothesis plays important role. This approach is 
alternative to the conventional one based on the dispersion 
relation for the Q^ dependent amplitude axial current in- 
teraction. While the latter faces the problem of lacking 
experimental information on most of the diffractive diago- 
nal and off-diagonal amplitudes, the dipole formalism is free 
of these difficulties. Besides, one can employ the universal 
dipole cross section (see Eq. ([3])), well fixed by numerous 
data for interactions of the vector current in electromag- 
netic processes (DIS, photoproduction, etc.). 

The important challenge of the dipole description is the 
construction of the current distribution amplitudes at small 
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Appendix A: Distribution amplitudes in the instanton 
vacuum model 



In this section we present the results of calculation of 
the DAs for the axial current and pion performed in the 
instanton vacuum model (IVM). 
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FIG. 8: Left: Cross-section of diffractive neutrino-production da/dW as function of W , for fixed neutrino energies E,^. Right: the 
same cross-section (da/dW) weighted with the neutrino spectrum from Minerva [J], see Eq. 1)35^ for exact definition. 



1. Axial DAs 

There are eight independent axial DAs defined in (jlSHlSp , 
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FIG. 9: Comparison of the color dipole prediction with experi- 
mental data from \5^. "CD" stands for the color dipole model, 
either with or without reggeons, "Adler relation" stands for 
the evaluation using Adler relation fl]) extrapolated to nonzero 
Q^ using dipole-type dependence for the Q^-dependence, ~ 



After tedious but straightforward calculations we arrive at. 
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2. Pion DAs 

For the pion there are four independent pion DAs defined 

in (EBlMl), 
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Eventually we arrive at the following structures in the pion 
DA, 
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The details of calculation of the distribution amplitudes 
will be presented elsewhere. 



3. Vector current DAs 



The vector current DAs were derived in |64| . 
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where the distribution amplitudes (^||,(^±7 have twist 2, The formfactors /^(g), /" (g), /^ (g^) and normahzation 
V'±^)V'7°\V'7'';2?T has twist 3; and /i^"-*, /i^*^ have twist 4. constants f3^,Xm,fA are discussed in detail in [63|. 
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